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ABSTRACT 

We have investigated the ionisation structure of the post-shock regions of magnetic 
cataclysmic variables using an analytic density and temperature structure model in 
which effects due to bremsstrahlung and cyclotron cooling are considered. We find that 
in the majority of the shock-heated region where H- and He-like lines of the heavy 
elements are emitted, the collisional-ionisation and corona-condition approximations 
are justified. We have calculated the line emissivity and ionisation profiles for Iron 
as a function of height within the post-shock flow. For low-mass white dwarfs, line 
emission takes place near the shock. For high-mass white dwarfs, most of the line 
emission takes place in regions well below the shock and hence it is less sensitive to 
the shock temperature. Thus, the line ratios are useful to determine the white-dwarf 
masses for the low-mass white dwarfs, but the method is less reliable when the white 
dwarfs are massive. Line spectra can, however, be used to map the hydrodynamic 
structure of the post-shock accretion flow. 
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1 INTRODUCTION 

When a supersonic accretion flow abruptly becomes sub- 
sonic near the surface of a compact star, a stand-oiT shock is 
formed. For the white dwarfs in magnetic cataclysmic vari- 
ables (mCVs) (see Warner 1995), the accretion-shock tem- 
perature Ts = 3 Gfj,m^Mv, /SkRv, (if Xa << Rw), where G is 
the gravitational constant, k the Boltzmann constant, 
the hydrogen mass, fi the mean molecular weight of the gas, 
Xa the shock height, Ad„ the mass of the white dwarf and _Rw 
the radius of the white dwarf. (Typically, Ta ~ 10® K.) As Ta 
depends on Mw and R^,, and is insensitive to other orbital 
parameters, the white-dwarf mass can be determined if Ta 
is measured (Rothschild et al. 1981; Kylafis & Lamb 1982, 
Ishida 1991; Fujimoto & Ishida 1997; Cropper, Ramsay & 
Wu 1998). 

The distance from the shock front to the white-dwarf 
surface (i.e. the shock height) and the structure of the shock- 
heated emission region depends not only by Ts but also the 
cooling processes. For mCVs, bremsstrahlung and cyclotron 
cooling are the most important. When the accretion rate of 
the system is high {i^ 10^® g s~^) and the white-dwarf mag- 
netic field is weak (~ 1 MG), bremsstrahlung cooling usually 
dominates. If the accretion rate is low (^ 10^^ g s~^) and 
the magnetic field is strong (^ 10 MG), then cyclotron cool- 



ing is the dominant process, at least near the shock (Lamb 
& Masters 1979; King & Lasota 1979). Line cooling is gen- 
erally unimportant, except at the bottom of the post-shock 
region, where the electron temperature Tc falls well below 
10^ K (e.g. Kato 1976; Mewe, Gronenschild & van den Oord 
1985). 

It is straightforward to calculate the post-shock accre- 
tion fiow with only optically thin bremsstrahlung cooling 
(see e.g. Aizu 1973; Chevalier & Imamura 1982). The inclu- 
sion of cyclotron cooling complicates the calculations be- 
cause of large opacity effects. An exact treatment of cy- 
clotron cooling requires the solving of the radiative-transfer 
and hydrodynamic equations simultaneously, and it is un- 
likely that simple analytic solutions can be obtained. How- 
ever, for parameters typical of mCVs, the total cooling pro- 
cess can be roughly approximated by a composite cool- 
ing function with power laws of density and temperature 
(Langer, Chanmugam & Shaviv 1982; Wu 1994; Wu, Chan- 
mugam & Shaviv 1994). With this composite cooling func- 
tion, the hydrodynamic equations governing the post-shock 
fiow can be solved, yielding a closed-form solution. With 
the density and temperature structures of the emission re- 
gion specified, the X-ray emission from the system can be 
calculated. (See Wu 2000, for a review of accretion shocks 
in mCVs.) 



© 0000 RAS 



2 K. Wu, M. Cropper & G. Ramsay 



X-rays emitted from regions in which both cyclotron 
and bremsstrahlung cooling are important are generally 
softer than those from regions with bremsstrahlung cooling 
only. Moreover, the spectra of the cyclotron-bremsstrahlung 
cooling shocks arc richer in emission lines (Cropper, Ramsay 
& Wu 1998; Tennant et al. 1998). Emission lines can also be 
used to diagnose the density and temperature structures of 
the shock-heated region (sec Fujimoto & Ishida 1995; Ezuka 
& Ishida 1999). The grating instruments on board the cur- 
rent X-ray satellites XMM-Newton and Chandra can easily 
resolve many of the emission lines in the ~ 1 keV range. 
With the spectral resolution of these satellites we will be 
able to measure the accretion flow very close to the white- 
dwarf surface. 

From the X-ray continuum, the shock temperature can 
be determined, and hence the white-dwarf mass (e.g. Ishida 
1991; Cropper, Ramsay & Wu 1998; Ramsay 2000; Beard- 
more, Osborne & HcUior 2000). In analysing X-ray data, 
model continua spectra arc fitted to the data in order to 
determine the system parameters. As noted above, numer- 
ical calculations with exact treatment of cyclotron cooling 
will in principle produce more accurate model spectra, but 
its heavy demand on computation time makes it imprac- 
tical in most situations. Semi-analytic methods that can 
yield reasonably accurate results are therefore more useful 
in the data analysis to extract system parameters, such as 
the shock temperature. In this study we investigate the ion- 
isation structures and the line emission of the post-shock 
emission regions of mCVs in terms of the analytic treat- 
ment of Wu (1994) and Wu, Chanmugam & Shaviv (1994). 
We ignore the effects of gravity over the height of the shock 
(e.g. Cropper et al. 1999) and unequal ion and electron tem- 
peratures (e.g. Imamura et al. 1987; Saxton 1999; Saxton & 
Wu 1999). 



2 STRUCTURED SHOCK-HEATED REGIONS 

The ionisation of the elements in the post-shock flow in 
mCVs is generally caused by X-ray irradiation or collisions. 
The electron transitions between the K-, L-shells and the 
outer shells are the major processes in producing the keV 
emission lines. The strength of the lines are therefore de- 
termined by the temperature, density, and the ionisation 
structures of the emission region. 

The electron temperature at the accretion shock is ~ 
10* K. This temperature is sufficient to completely ionise 
the elements such as Argon, Silicon, Sulphur, Aluminium or 
Calcium through electron collisions. Iron can be ionised to 
the highest ionisation states: Fe XXVI and Fe XXVII. 

The ionisation due to irradiation can be described by 
the ionisation parameter, which is defined as ^ = Ly^/n^r^ 
(Kallman & McCray 1982). (Hero, Lx is the luminosity of 
the irradiation, is the rmmbcr density of hydrogen, and r 
is the characteristic size of the plasma.) In the shock-heated 
region of mCVs, Lx ~ 10^^ - 10^^ erg s'^ , ~ lO^'^ - 
10^® cm~^, and r ^ 10* cm and the ionisation parameter is 
^ lO'^ . The corresponding ionisation states of Iron are Fe XX 
or lower (Makishima 1986). 

Therefore, in the shock-heated regions of mCVs colli- 
sional ionisation is the major process producing the keV 



emission lines. Photo-ionisation is important only in the 
cooler and less dense pre-shock accretion stream. 

2.1 Collisional ionisation equilibrium 

For T ~ 10^ — 10* K, collisional ionisation equilibrium can 
be attained when net > 10^^ cm~^s, where Ue is the electron 
number density and t is the time spent by the ions in the 
plasma (Gorenstein, Harnden & Tucker 1974; Masai 1984). 
For a white dwarf with a specific accretion rate m, the local 
value of riet at C, (the distance normalised to the shock height 
Xs) is 



m f 



v' dv' 



1 



2GM^m„ T 



1/4 



dr' 



dr' 



(1) 



where t = —v/v^^ is the dimensionless velocity, and % is 
the free-fall velocity at the white-dwaxf surface. 

To evaluate the integral in the above equation we con- 
sider the inverted velocity profile of the accretion fiow given 
in the liydrodynamic model of Wu (1994) (see Cropper, 
Ramsay & Wu 1998 for typographical corrections): 



dr 



r^(5-8r) 



1 + 3-"4"+''' 



(2) 



where K{t) = 1 + 3-"4"+'^es(l - r)°r'' (with a « 2 and 
P ~ 3.85 appropriate for cyclotron cooling), ps is the density 
at the shock surface, ts is the ratio of the bremsstrahlung- 
cooling time-scale to the cyclotron-cooling time-scale at the 
shock, and A — 3.9 x 10^^ in c.g.s. units. The shock-height 
Xs can be obtained by a direct integration of the inverted 
velocity profile, i.e., 

»l/4 



Aps 



dr- 



'(5-8r) 



K{r) 



(3) 



In this formulation, the hydrodynamic variables, such as 
density p, pressure P and temperature T, are proportional 
to 1/r, (1 — r), and r(l — r) respectively (see Wu 1994; Wu, 
Chanmugam & Shaviv 1994). 

By combining Equations (1) and (3), this yields 

[net] w2.2 X 10^^ Fi(C;£s) 



Mo 



1/2 



5 X 10* cm 



-1/2 



-3 

cm s 



where 



■8r') 



Kir') 



(4) 



(5) 



As t ^ Xs/v^f ^ min(l,e7^), when eg is specified [nct]^ 
does not explicitly depend on r'n (sec Equation 4). In Fig- 
ure 1, wo show 2.2 X lO** -Fi(C; £s) as a function of (1 — (^), the 
distance from the shock surface. When 2.2 x 10* i^i(C; Eb) > 
1, the accreting matter is in collisional ionisation equilib- 
rium. Figure 1 shows that except in a very thin layer near 
the shock front itself, collisional ionisation equilibrium is 
reached in the post-shock region for the range of eg of inter- 
est. 
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Figure 1. The quantity 2.2 X 10'* Fi (C; Cs) as a function of (1 - C), 
the normahsed distance from the shock surface. The plasma is in 
coUisional equilibrium if 2.2 X 10*Fi(^; eg) > 1. Curves a, b, c and 
d correspond to es = 0, 1, 10 and 100 respectively. 
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Figure 2. 5 x 10^ Fi(C,; es, Z) as a function of C, the normalised 
height from the white-dwarf surface, for Iron (Z = 26). For the 
solid curves Mw = I.OMq and = 5.53 x 10* cm; and for 
the dotted curves, Mw = O.5IVI0 and R^r = 9.61 x 10* cm. In 
each set of the curves, the corresponding are 0, 1, 10 and 
100 (from top to bottom respectively). The horizontal line repre- 
sents F2{(^; ee, Z) = 2.0 x 10"^, above which the approximation 
of corona condition is applicable. 



2.2 lonisation balance 



The lonisation balance is determined by the excitation and 
de-excitation rates of the corresponding lonisation states. 
As the rate coefficients may depend on the optical depth, 
it is generally required to solve the coupled population rate 
and radiative-transfer equations to obtain the population of 
the lonisation states (Bates, Kingston & McWhirter 1962; 
Castor 1993). However, if the electron number density is 
sufficiently low that the radiative-transfer effects arc unim- 
portant, the lonisation balance can be calculated using the 
corona-condition approximation. In the other extreme, if the 
electron number density is sufficiently high such that the 
plasma is able to attain local thermal equilibrium (LTE), 
the population can then be determined by the Saha equa- 
tion. 

For the corona condition to hold, it requires the gas 
to be in coUisional equilibrium and a direct coupling of the 
heat input to the ions and the free electrons. It also requires 
the density of the plasma to be sufficiently low such that 
most ions are in ground states, and the line (and continuum) 
emission to be optically thin (Elwert 1952; Mewe 1990; Kahn 
& Liedahl 1991). 

For ions with a bared nuclear charge Z, the corona con- 
dition is satisfied if rie ^ 4 x IQ^Z^ cm"* (Wilson 1962; 
Mewe 1990). As rie and Te scale with 1/r and t(1 — r), 



where r is the dimensionless velocity (see §2.1), we can de- 
fine a quantity 



F2(C;es;2)=^MCr[l-r(C)]' 77^ 



Mq 



5/2 



0.5 / V 1 g s 



-5/2 



5 x 10* cm 



(6) 



The corona-condition criterion is then equivalent to require 
F2iC,e.;Z) > 2 x 10-^ 

In Figure 2 we show the profile of the function [5 x 
lO^F^iC; Es, Z)] in the shock-heated regions of accreting mag- 
netic white dwarfs for various parameters. It can be seen that 
for white dwarfs with Mw ~ 1 M© and m = 1 g s~^, the 
corona condition is satisfied when the height ^ > 10~^. For 
white dwarfs with Mw = 0.5 M©, C > 10"*. The corona- 
condition approximation is not applicable at the very bot- 
tom of the shock-heated region, where the electron number 
density is high and the temperature is low. 

A plasma is in LTE when ne ^ 1.4 x lO^^Z^T^/^ cm"* 
(Wilson 1962). In terms of the velocity profile of the post- 
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shock accretion flow and white-dwarf parameters, it requires 
1 > 9.7 X 10^ZV(C)^/^[1 - r(C)]^/' {j^ 

^ (l g s-i) (m^) (5 X 108 cm) ■ ^'^^ 

For Iron, Z = 26. LTE implies t ^ 10~^, corresponding to 
t; ^ 1 cm s^^. Thus, LTE is applicable for the very base of 
the post-shock flow only. 

2.3 lonisation profiles 

lonisation-balance calculations in the approximation of 
corona condition have been carried out by many authors 
(e.g. Jordan 1969; Jacobs et al. 1977; Raymond & Smith 
1977; Mcwc & Groncnschild 1981; Arnaud & Rothcnflug 
1985; Mewe, Gronenschild & van den Oord 1985). The re- 
sults of the calculations are usually presented in terms of 
riz/ riz (the relative ion concentration) as a function of 
Te (the electron temperature). With the local value for the 
electron number density, the electron temperature and the 
abundance of the elements specified, the results of these ion- 
isation balance calculations can be readily used to determine 
the concentration of the ion species. 

The temperature of the accretion shock in mCVs is typi- 
cally ~ 10® K, which is sufficient to ionise Iron in the shock- 
heated region to its higher ionisation states. As only the 
most populous ion species would make significant contribu- 
tion to the emission, in this study wo only consider the four 
most highly ionisation states, Fe XXVII {¥c^'^+), Fe XXVI 
(Fe25+), Fe XXV {¥e'^'^+) and Fe XXIV (Fc2-^+). Also, we re- 
strict the emission region of interest to be C > 10"'', where 
the corona condition is satisfied. 

Wc first determine the local ionisation by interpolating 
the results from the ionisation balance calculations of Ar- 
naud & Rothenflug (1985) and then use it to calculate the 
profiles of the relative concentration of the ion species in 
the shock-heated region. Two representative cases are con- 
sidered, with the white-dwarf masses 1.0 and 0.5 re- 
spectively for each case. (We have Eissumed the Nauenberg 
(1972) mass-radius in the calculations.) The parameter £a is 
chosen to be 0, 1, 10 and 100. This covers the range of pa- 
rameters appropriate for most mCVs (see Wu, Chanmugam 
& Shaviv 1994). 

In Figures 3 and 4 we show the relative ion concentra- 
tions of Iron as a function of the height above the white- 
dwarf surface for the I.O-M0 and O.5-M0 white dwarfs re- 
spectively. The normalisation is such that the total concen- 
tration is 1.0 at the shock surface {C, = 1). As the density 
(and the electron number density rie) increases when ^ de- 
creases, the peak concentration of each ion species is always 
> 1.0. Since the temperature and density in the post-shock 
region arc monotonic functions of C, (if the electrons and ions 
have the same temperature and the variation of gravity with 
height is negligible), when the temperature is specified, the 
density is uniquely determined. In ionisation equilibrium, 
the peak concentration of the ion species depends only on 
Te. Thus, for a fixed white-dwarf mass, the peak concentra- 
tions occur at the same Te and rie, in spite of the different 
values of e^. 
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Figure 3. The ionisation profile of Fc XXVII, Fc XXVI, Fe XXV 
and Fc XXIV in the sliock-licatcd emission region of a I.O-M0 
wfiite dwarf. Tfic effects of cyclotron cooling increases from the 
top to the bottom panels (from panel (a) to panel (d) respec- 
tively), and the corresponding parameters are €3 = 0, 1, 10 and 
100. The surface of the white dwarf is at C = 0; and the shock 
surface is at C = 1- 
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Figure 4. Same as Fig. 3 for a white dwarf with 0.5 Mq. 
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Figure 5. he mean charge profile of Iron in the shoclc-heated 
region for (a) Mw = I.OMq; and (b) for Mw = O.SMq. Curves a, 
b, c and d correspond to eg = 0, 1, 10 and 100 respectively. 

As shown in Fig 3, Fe XXVII dominates in most of the 
shock-hcatcd region of the I.O-M0 white dwarf with eg = 0. 
The lower ionisatiori states are important only near the low- 
tcrriperature bottom of the post-shock region. For larger es, 
the emission region is cooler, and so the lower ionised species 
becomes more abundant. When es = 100, the relative con- 
centration of Fe XXVI is higher than that of Fe XXVII in the 
region where ~ 0.1 — 0.3. However, as a whole, Fe XXVII 
still dominates in the post-shock region. For a 0.5-M© white 
dwarf, the shock temperature kTg fs 13.6 keV, which is 
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Fe XXVI [HI] 
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0.01 

normollsed height ^ 

Figure 6. The emissivity (in erg cm~^ profile of the 

Fe XXVI Lym£ui-a (Is -^S - 2p ^P) transition for a white dwarf 
with a mass of 1.0 Mq (upper panel) and 0.5 Mq (lower panel). 
The specific accretion rate is m = 1 g cm~^ s~^. Curves a, b, c 
and d correspond to eg = 0, 1, 10 and 100 respectively. For com- 
parison, wc also show in dotted line the line emissivity of emission 
regions with constant electron temperature and number density 
the same as the values at the shock of the structured emission 
regions that we consider. 



insufficient to completely ionise all the Iron to Fe XXVII. 
For es^l the dominant species is Fe XXVI, and for es^W, 
Fe XXV. As in the previous case, the lower ionisation states 
becomes more important at the bottom of the post-shock 
region. 

In Figure 5, we show the mean charge (z) profiles. For 
the I.O-M0 white dwarfs, z ~ 25.8 at the shock surface, 
which implies that Fe XXVII is the most abundant ion. As 
C decreases, z first decreases slowly to « 24 and then falls 
more rapidly afterwards. For es = (bremsstrahlung cooling 
only), z ~ 23.2 at = 10"''; while for es — 100 (cyclotron 
cooling dominated), z ~ 18 at = 10"'*. For the 0.5-Mq 
white dwarfs, z « 25.3 (i.e. Fe XXVI dominates) at C = 1-0. 
Similarly to the l.O-M© case, z decreases slowly with ^ until 
it reaches « 24 and then rapidly afterwards. For es = 0, 
« = 17 at C ~ 10"*; and for = 100, « = 17 at C ~ 10"^. 
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Figure 7. Same as Fig. 6 for the Fe XXV He4 (Is^ - ls2p 
^P) transition. 



2.4 Line emissivity profile 

In Figures 6, 7, 8 and 9, the line emissivity of the Fe XXVI 
Lyman-a (Is - 2p ^P), Fe XXV Hc4 (Is^ 'S - ls2p ^P), Fe 
XXV He5 (Is^ - ls2p ^P), and Fe XXV He6 (Is^ - ls2p 
^S) lines are shown. The corresponding line centre energies 
of these lines are 6.93, 6.70, 6.68 and 6.64 keV. In calculating 
the emissivity we have interpolated the table of line power 
given in Mewe, Gronenschild & van den Oord (1985). The 
electron number density that we assume corresponds to an 
specific accretion rate of m = 1 g cm"^ s"^. 

'When the white-dwarf mass and the accretion rate are 
fixed, the emissivity of an emission line is the same at the 
shock for all eg. As ts increases, cyclotron cooling becomes 
more efficient. As the electron temperature of the shock- 
heated region is lowered, the lower ionised species becomes 
more significant. As a result, the cmissivities of the Fe XXVI 
and Fe XXV lines all peak at larger ^ for larger ts (Fig. 10). 



3 DISCUSSION 

3.1 White-dwarf mass determination 

"When the local emissivity Pi(C; Ts, Cs, z\) of an emission line 
z\ is specified, the intensity of the line is simply the sum of 
contribution from all zones in the emission region, i.e., 
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Figure 8. Same as Fig. 6 for the Fe XXV He5 (Is^ - ls2p 
^P) transition. 
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Figure 9. Same as Fig. 6 for the Fe XXV He6 (Is^ - ls2p ^S) 
transition. 



x^S 



A{Z) / dCPi{C;T,,e.,zO 



(8) 



where A{Z) is the abundance of the element that give rise 
to the emission hne Zi, Ai the height above the white-dwarf 
surface at which the emission becomes optically thick, S 
the cross-section area of the emission region, and D is the 
distance of the source. In terms of the temperature gradient 
dT/dC, 



Ta, 



dT^n{C;T,,e,,zO , (9) 



where Taj is the temperature at Ai. The ratio of the inten- 
sities of two emission lines Zi and Zj of the same element Z 
are therefore 

ItI dT§m;Ts,e.,zO 

Rizi,zy,Ts,e,,Ai,Ai) = — -. (10) 

J^l^dT§Pi{C;Ts,e,,Zi) 

As the line-intensity ratio is a function of ts, 2s, Taj and 
Taj only, one may use the observed line intensity ratio to 
constrain these parameters. 

Fujimoto (1996) and Fujimoto & Ishida (1995, 1997) 
assumed that Taj = 2Ai(= Ta), and use the ratio of El- 
and He-like Fe Kq lines observed by ASCA to determine 
the shock Tg of the Intermediate Polar EX Hya. Prom the 



deduced shock temperature they obtained a white-dwarf 
mass of 0.48lg jg M0. Using the same method, Hellier et 
al. (1996) deduced the white-dwarf mass of another Inter- 
mediate Polar AO Psc, giving 0.40l|',;"^ Mq. These masses 
arc generally in agreement with those determined by fit- 
ting the X-ray continuum using Ginga data (0.45 Mq and 
0.45 Mq respectively for EX Hya and AO Psc) and RXTE 
data (0.44±0.03 Mq and 0.60±0.03 Mq respectively) (Crop- 
per, Wu & Ramsay 1999; Ramsay 2000). 

As the emissivity of the emission lines is sensitive to the 
temperature of the emission region, it is generally considered 
that white-dwarf masses can be accurately determined us- 
ing spectral lines, and the line emission may provide better 
constraints to the white-dwarf masses than the continuum. 
However, an accurate measurement of the line strengths does 
not necessarily lead to an accurate determination of the 
shock temperature. The shock temperature can be deduced 
only if there are sufBcient concentrations of ions that are re- 
sponsible for the line omission at the shock. If the shock tem- 
perature is sufficiently high that the element is completely 
ionised at the shock, the strengths of its H- and He-like 
lines are insensitive to the shock temperature. These lines 
arc now omitted from the cooler bottom of the post-shock 
region (cf. Fig. 10a and Fig. 10b). Thus, using Fe H- and He- 
like emission lines to determine white-dwarf mass is practical 
for the low- mass (~ 0.5 Mq) systems but not for massive 
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Figure 10. (a) The "(^ x Pi" emissivity profile of the continuum 
at 6.8 keV (dotted curve), the Fe XXV He4 hne (curve a), the 
Fe XXV He5 hne (curve b), the Fe XXV He6 line (curve c) and 
the Fe XXVI Lyman-a line (curve d) for a 0.5-Mq strongly mag- 
netised white dwarf with es = 100. (b) Same as (a) for a l.O-M© 
white dwarf with €« = 0. 



white dwarfs (Wu 2000). Indeed, an attempt to constrain 
the white-dwarf masses of two more massive AM Her type 
systems, AM Her itself and BL Hyi, using this method was 
unsuccessful (Fujimoto & Ishida 1995). Although one may 
consider elements with higher iouisation potential, such as 
Nickel, the abundance of heavier elements in typical mCVs 
are generally low and unable to produce lines that are strong 
enough to be useful in typical observations. For instance in 
Figures 11, 12, 13 and 14 we show the strengths of the Fe 
and other lines in simulated XMM-Newton EPIC-PN and 
RGS spectra of accreting white dwarfs of different masses 
and cooling efficiencies. In these simulations we have used 
calibration files determined 7 months into flight and used 
exposures typical of XMM-Newton observations. Tennant et 
al. (1998) shows simulated spectra using the Chandra grat- 
ing spectra. 



3.2 Flow diagnosis 

Because of the temperature and velocity stratification in the 
shock-heated region, the emission lines from different heights 
above the white-dwarf surface have diflferent Doppler shifts. 




1 2 
chonnel energy (keV) 



XMM-Newton EPIC-PN: 100 kseo 



> * — 
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chonnel energy (keV) 



10 



Figure 11. (a) Simulated XMM-Newton EPIC-PN spectra of a 
0.5-Mq white dwarf with substantial cyclotron cooling (es = 10) 
for a 100-ks exposure. We assume solar metal abundance. The 
flux is normalised to that of the observed flux of AM Her itself 
in the 2— lOkeV energy band (fa 7 X 10~^^ erg s~^ cm~^) derived 
from ASCA data (0.5-lOkeV) (Ishida et al. 1997). (b) Same as (a) 
for a I.O-Mq white dwarf with no cyclotron cooling. In simulating 
the spectra we have used calibration files derived 7 months into 
the orbit of XMM-Newton and assume a medium filter. 



A line can be emitted only if the temperature of the ac- 
cretion matter allows the according electronic transition to 
occur, so one can infer the temperature and density of the re- 
gion where the line is emitted. In addition the Doppler shift 
of the line centre energy can be mecisured, from which the 
line-of-sight velocity of the emitter can be deduced. Thus, 
by examining the different lines in a spectrum, one can ob- 
tain a relation between the flow velocity, temperature and 
density. 

For a l.O-M© white dwarf, the velocity of material at 
the accretion shock is about 2 x 10^ cm s~^. The centre en- 
ergy of lines emitted from regions near the shock will have 
an energy shift of about 0.6%. As the white dwarf in a mCV 
revolves with the orbit and it also rotates, there is an addi- 
tional velocity shift of about 2 — 3 x 10^ cm s~^ for typical 
CV parameters. The white-dwarf rotational period and the 
binary orbital period are the same for typical AM Her type 
systems, and it allows us to subtract the velocity due to or- 
bital motion and the white-dwarf rotation easily. Moreover, 
AM Her type systems do not have an accretion disc, and 
therefore the line spectra are emitted mainly from the hot 
shock-heated accretion matter near the white-dwarf surface. 
The Doppler shift of the emission lines is then only due to 
the bulk accretion flow, and they can be extracted as de- 
scribed above from the orbital phase-binned data. 
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Figure 12. (a) Simulated XMM-Newton RGS spectra (0.3- 
2.1keV) of a 0.5-Mq wliite dwarf witli substantial eyelotroii cool- 
ing (cs = aO) for a 100-ks exposure. We assume solar metal abun- 
dance. The flux is normalised to the observed flux of the hard 
component of AM Her itself (» 1.7 x 10~^^ erg cm~^) deter- 
mined from i?05AT data (0.1-2.4keV) (Ramsay et al. 1994). (b) 
Same as (a) for a l.O-M© white dwarf with no cyclotron cooling. 
In simulating the spectra we have used calibration files derived 7 
months into the orbit of XMM-Newton. 



A liric-ccntrc energy resolution of 1 part in 10000 
(achievable by XMM-Newton and Chandra at about 1 keV) 
will allow us to measure the accretion flow a few tens metres 
above the white-dwarf surface. As shown in Figures 13 and 
14, the Fc complex and other lines arc clearly resolved in 
the XMM-Newton simulated grating spectra for a 100-kscc 
observation. In the parallel study of Tennant et al. (1998), 
it was shown that an energy resolution of 0.1-eV can be 
achieved by the Medium Energy Gratings of Chandra for a 
single gaussian fit to a line of centre energy of 1 keV for 
about 200 photons in the linos. For a mCV of the X-ray 
brightness of the system AM Her itself, this requires an ob- 
servation of about 100 ksec to detect this number photons 
in 1/10 of an orbital phase. 

X-ray spectroscopy will therefore allow the flow velocity 
in these accreting systems to be measured directly. This di- 
agnostic power cannot be achieved by optical spectroscopy. 



despite its better velocity resolution, as optical emission 
from the shock-heated material in mCVs is optically thick. 



4 SUMMARY 

We have investigated the ionisation structure and line crnis- 
sivity profiles of the shock-heated emission regions of accret- 
ing white dwarfs in mCVs by means of an analytic model 
given in Wu (1994) and Wu, Chanmugam & Shaviv (1994). 
We have found that for Iron, the corona-condition approx- 
imation is generally satisfied in most of the shock-heated 
region where lines are emitted. Because of the temperature 
and density stratification in the shock-heated region, differ- 
ent lines arc emitted from different heights above the white- 
dwarf surface. By measuring the Doppler shifts of the lines, 
one can obtain a relationship between the flow velocity, the 
temperature and the density of the emitting gas, thus pro- 
viding a means to map the hydrodynamics of the post-shock 
accretion flow directly. Our study also indicates that using 
emission lines to determine white-dwarf masses is practical 
only for low-mass systems. For massive systems, the abun- 
dant elements such as Fc are completely ionised at the shock 
and so their line emission is insensitive to the shock temper- 
ature. 
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